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We study, by means of the stabilizer formalism, a quantum error correcting code which is alter- 
native to the standard block codes since it embeds a qubit into a qudit. The code exploits the 
non-commutative geometry of discrete phase space to protect the qubit against both amplitude 
and phase errors. The performance of such code is evaluated on Weyl channels by means of the 
entanglement fidelity as function of the error probability. A comparison with standard block codes, 
like five and seven qubit stabilizer codes, shows its superiority. 

PACS numbers: quantum error correction (03.67.Pp); decoherence (03.65. Yz) 

I. INTRODUCTION 

In order to describe realistic quantum information processes, quantum errors induced by environmental noise must 
be taken into account. This can be accomplished by introducing the notion of quantum channels, that is to say maps 
on the set of states of the system that are completely positive and trace preserving [T]. At the same time, one would 
combat quantum errors to avoid their detrimental effect on quantum information processes. To this end, the method 
of error correcting codes has been borrowed from classical information theory (for a comprehensive introduction to 
the quantum theory of error correcting codes we refer to [2]). The underlying idea is to exploit redundancy, that is to 
encode information in linear subspaces (codes) of the total complex Hilbert space in such a way that errors induced 
by the interaction with the environment can be detected and corrected. 

Usually a logical qubit (a two dimensional complex Hilbert space H2) is encoded into n physical qubits (a 2™ 
dimensional complex Hilbert space Wf ™). This kind of encoding is known as block- encoding 

7~l~2 ^ Illogical) 1 > | (/physical) G Wf™. 

However, there is also the possibility of embedding a logical qubit into a d-dimcnsional quantum physical system, i.e. 
a qudit with complex Hilbert space T-Ld ^ Hf 71 - We refer to this kind of encoding as embedding or qudit- encoding 

n 2 9 I ^logical ) 1 > ^physical) € %d ^ %2 ™, 

where d ^ 2™ . 

For block-coding schemes a powerful formalism, named stabilizer formalism [3], has been developed describing one 
of the most important classes of quantum codes, namely the quantum version of linear codes in classical coding theory. 
The stabilizer formalism can be extended over non binary codes [J] and can be also useful for describing embedding 
(qudit) codes. Actually, the idea of embedding a qubit into a larger space without resorting to block codes was put 
forward in [5] where a qubit was encoded into a bosonic mode (infinite dimensional complex Hilbert space) just using 
the stabilizer formalism. Later on the possibility of qudit encoding was pointed out in Ref. [B] by using the same 
formalism. There, being in "Hd, the errors were considered as a generalization of Pauli operators, representing shift 
errors (A-type) or phase errors (Z-type) or their combination (AZ-type). However, in [6] the proposed code was 
essentially classical since only Z-type errors were taken into consideration. 

In this article we upgrade such a coding scheme to be fully quantum, thus able to correct A-type, Z-type and, 
AZ-type errors. We then test its effectiveness for d = 18 and d = 50 on a Weyl quantum noisy channel [7] (an error 
model characterized by errors of A, Z and AZ-types). We allow for the possibility of considering A and Z errors 
occurring with both symmetric and asymmetric probabilities. Finally, we compare the performance of such qudit 
coding schemes to those of the conventional five-qubit stabilizer code [[5, 1, 3]] [8, 9] and seven-qubit stabilizer CSS 
(Calderbank-Shor-Steane) code [[7, 1, 3]] [TUHH]. We characterize the performances of these codes by means of the 
entanglement fidelity |12) , rather than the averaged input-output fidelity used in [6] . We show that the use of qudit 
codes may allow to save space resources while achieving the same performance of block codes. 

The layout of the article is as follows. In Section II, we briefly describe block-encoding and then introduce qudit 
encoding. Special focus is devoted to both Pauli groups of n-qubit vectors and generalized Pauli groups of error oper- 
ators acting on a qudit. In Section III, the Weyl noisy quantum channel is discussed together with the entanglement 
fidelity. In Section IV, we study the performance of qudit codes. In Section V, for the sake of comparison, we quantify 
the performance of relevant block codes. Our final remarks appear in Section VI. 
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II. FROM BLOCK CODES TO EMBEDDING CODES 

In this Section, we briefly recall the block-encoding error correction schemes in terms of the stabilizer formalism. 
Then, using the same formalism we introduce the qudit-encoding (embedding qubits into qudits) scheme. In both 
cases we restrict our attention to a single encoded qubit. 



A. Pauli group of n-qubits and block-encoding 



A qubit is a two-dimensional quantum system with associated complex Hilbert space H.2- Let {|0),|1)} be the 

def t r def . r? def i n i • r 

Y = ic&y, Z — c 2 denned m terms ot 



denotes the imaginary 



canonical basis of this space and consider on it the Pauli operators X = a 
the standard Pauli operators o~ x , o~ y , a z realizing the su(2) algebra (throughout the paper 
unit of C). They are such that {|0), |1)} are eigenstates of Z, 

xz = -zx, 

and Y = XZ. The Pauli operators so defined suffice to describe all possible errors occurring on a single qubit. 
Together with the identity operator / (with / = I2X2) they form a multiplicative group if we allow them to be 
multiplied by —1, i.e. {±J, ±X, ±Y, ±Z}. We refer to this group as the Pauli group Vn 2 - Actually it is a subgroup 
of the Pauli group realized through the standard sigma Pauli operators and it coincides with the discrete version of 
the Heisenberg-Weyl group. 

For n-qubit errors we can then consider the Pauli group Vn,®n whose elements result from ra-fold direct products 
(see also [T3]) 

e(A,ji, . . . , j n ) = (-) A e 1 (j 1 ) ® ... ® e„(j„). 

The subscripts on the RHS label the qubits 1,..., n, while jk = 0,1,2,3, label respectively the operators Ik, Xk, Yk, 
Zk acting on the fc-th qubit. Furthermore A £ {0, 1}. 

Since efc(2) = XkZ^, the elements e(A, ji, . . . , j n ) can be rewritten as 

e(\,a,b) = (-) x X(a)Z(b), 

where a = a\...a n and b — bi...b n are bit strings of length n and, 



X (a) d ^ f (A!) ai <8 - ® {X n ) an , Z (b) a ^ {Z,p ® ... ® (Z n ) 



def 



Observe that the order of the group P M «n of n-qubits errors is = 2 2n+1 . Since the factor ± in front of an 

error makes no relevant difference on its action, we can actually assume to work with the quotient group V H ® n / {±/}, 
with the major exception being the determination whether elements commute or anti-commute. 



The order of the quotient group is 



7W{±I} 



y2n 



Then, there is a one-to-one correspondence between 

■P w ®n/ {±1} and the 2n-dimensional binary vector space F|" whose elements are bit strings of length In [H]. A 

vector v £ Fj™ is denoted v = (a\b), where a = a\...a n and b — bi„.b n are bit strings of length n. Scalars take 
values in the Galois field F2 = {0, 1} and vector addition adds components modulo 2. In short, we have the following 
correspondence e(A, a, b) £ V u ®-n. <h> v e = (a\b) £ F|™. 

A quantum stabilizer code C is a vector space C C "H^" stabilized by an Abelian subgroup S C V^<»n , i.e. such that 



A quantum stabilizer code C with stabilizer generators 



fJ 



(n-l) 



is a 2-dimensional code space, i.e. a space 



where to encode a single qubit. The codewords |0) l, (basis vectors for such code space) can be found as orthogonal 
eigenvectors (corresponding to the eigenvalue +1) of any of the generators j®, 
The encoding operation then reads 



0n 



H 2 3 |0) 1— > |0 L ) e Hf n and, H 2 3 |1) h> \l L ) € H 

In view of the correspondence between the Pauli group and the vector space F^™, let 
of the generators gW> in F 2 " and let introduce the so called parity check matrix 

( ( a m\bW) \ 
( a W\bW) 



H 



def 



(a^\b^) be the image 



(1) 



V (a^-^lb^-V) J 
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Then, for an error e € ■<-> w e = (a\b) € Ff", the error syndrome 5(e) is given by the bit string [T3] 

5(e) = ffu e = «i.../„_i, (2) 

where ^ = if T (j) • u e . 

Errors with non- vanishing error syndrome are detectable. They correspond to operators not in 5 and not commuting 
with those in 5. That is, a set of error operators £ C V u %n is detectable if £ ^ Z(S) — 5, with 2(5) the centralizer 

of the subgroup 5. Furthermore, the set of error operators £ C P^s- is correctable if the set given by £ ^ £ is in turn 

detectable [15], i-e. £ f £ ^ 2(5) - 5. 

It would be awfully tedious to identify either detectable errors or sets of correctable errors. However, the quantum 
stabilizer formalism allows to simplify such task [3]. This is a consequence of the fact that by means of such formalism 
it is sufficient to study the effect of the error operators on the generators of the stabilizer and not on the codewords 
themselves. Actually the syndrome extraction corresponds to measure the stabilizer generators. 

Finally, we denote by [[n, k, dc]] a quantum stabilizer code C with code parameters n (the length), k (the dimension) 
and dc (the distance) encoding k- logical qubits into rt-physical qubits and correcting [ rfc ~ 1 j-qubit errors ( [^J denotes 
the largest integer less than x). 



B. Generalized Pauli group and qudit-encoding 

A qudit is a c£-dimensional quantum system with associated complex Hilbcrt space Ud- On this space we can 
introduce a generalized version of the Pauli operators X and Z considered in the previous Subsection. They can be 
defined through their action on the canonical basis {|fc)}fcez d of %d [16] . 

X\k) = \k®l) and, Z\k) = uj k \k) , keZ d , (3) 

where "©" denotes addition of integers modulo d and w = exp (ic^f ) is a primitive dih root of unity (ui d = 1). The 
X and Z operators so defined are unitary (X* = X^ 1 and Z^ = Z^ 1 ), but not Hermitian, and satisfy X d = Z d = I 
(with / = Idxd) together with the commutation relations 

X a Z b = uj- ab Z b X a . (4) 

It is then possible to consider the Pauli group Vu d consisting of all operators e of the form 

e(l, n, m)=uj l X n Z m , 

where l,n,m € 7Ld- Similarly to the previous Subsection, for errors on a qudit we may refer to the quotient group 

Vu {--'II ') </ 1|- ' 

In addition to the reasons stated in the introductory Section about passing from block to embedding codes, another 
motivation is to understand whether or not finite dimensional versions of the shift-resistant quantum codes of Ref. [5 
are effective. Hence, following Ref. [SJ, we consider d = 2rir2, and introduce a code C stabilized by the Abelian 
subgroup 5 C Vu d generated by X 2r± , Z 2r2 . The codewords (basis for C) are eigenstates of Z 2r2 and X 2ri with 
eigenvalue one. Hence they only contain \k)s with values of k that are multiples of r\ and that are invariant under a 
shift by 2ri. They read 

|0L> = ^(|0) + |2r 1 ) + ... + |2(r 2 -l)r 1 )) , (5) 



' (|r 1 > + --- + |(2(r 2 -l) + l)r 1 » . (6) 



The encoding operation becomes in this case 

n 2 3 |0) k> |0 L > G U d and, H 2 3 |1) |U) e U d - 

If the states ([5|, (|6| undergo an amplitude shift, the value of k modulo T\ is determined by measuring the stabilizer 
generator Z 2r2 : and the shift can be corrected by adjusting k to the nearest multiple of r\. 

The codewords in the basis of X's eigenstates can be found by observing that the eigenstates of X and Z operators 
are connected by the Fourier transform 

w = ^X>-'>i;>, 

3=0 
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where X\i) = u l \i). Then, it turns out that 

|0 L ) = -L (|0) + K + • • • + |2(rpT)r 2 )) , (7) 

|U) = -L (R) + • • • + |(2(n ^TT+ l)r 2 )) . (8) 

The codewords Q, have the same form of §5§, ([6]), but with r\ and r 2 interchanged. Hence, if they undergo a 
phase shift, the value of k modulo r 2 is determined by measuring the stabilizer generator X 2ri , and the shift can be 
corrected by adjusting k to the nearest multiple of r 2 . 

To understand what is the set of correctable errors according to the condition £'£ ^ Z(S) — S we have with the 
help of Q 

f X a' z b'^ [X a Z h )X 2 ^ = e 2 ™( b - & ')A-2 X 2r, ( X a 'z b '\ {X a Z b ) , 

(x a ' Z b '^ (X a Z b )Z 2r2 = e -2^(a-a')M Z 2r2 (x a> ' Z h '\ (X a Z b ) . 

The phases on the right hand sides are non-trivial only if 

\a-a'\< 7 ±, and \b-b'\<^. (9) 

Therefore, the code C can correct all shifts corresponding to these conditions. They amount to r\r-2 — d/2 and this is 
also the number of possible error syndromes. 

Finally, notice that other families of qudit codes could be constructed by generalizing the Pauli operators in a 
different way, e.g. by making them Hermitian. 



III. THE WEYL QUANTUM CHANNEL AS ERROR MODEL 



In this Section, we first discuss the Weyl quantum channel for qudit states and for qubit states, and then the 
entanglement fidelity as quantifier of codes performances. 



A. General form of the Weyl quantum channel 

Consider a completely positive trace preserving map (CPT map or quantum channel) Aw, 

A< d > : 6 (H d ) A< d > (p) e & (H d ) , 

where & (Hd) is the set of positive unit trace linear operators in 7id- Then, A^ d > is called bistochastic if, 

A(d) fy«i. 



where / (= Idxd) is the identity operator in Hd- The d-dimensional Weyl channel is a bistochastic quantum channel 
of the following form [7], 



d-l 

A Weyl (P) = f 7r ( n ' m ) Un > ™P U n, mi ( 10 ) 



n, m=0 



where p £ & (Hd) and 7r(n, m) is an arbitrary probabil ity distribution, thus respecting < 7r(n, m) < 1, 
Sn _ m=o 7r ( n ' TO ) = ^ ne un itary Weyl operators U n> m in (10) are defined as, 



-l 



U n , m d = 5Z exp ( ic^jkm \ \k®n) 



dcf ■ 
fe=0 
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where "©", as specified earlier, denotes the addition of integers modulo d. They also satisfy the (Weyl) commutation 
relations Q, i.e. 



^n, rrJJ n' , m' exp 



■ 2?r I I >\ 

zc — (n m — nm ) 
a 



where < n, n! , m, m' < d — 1. Notice that U Uj m may be rewritten as U n> TO = J7„ ! oC^o, m- Furthermore, U n fi = X r 
and J7o,m = Z m . Therefore in what follows we will consider the d-dimensional Weyl channel acting as 



d-l 



A Wey» = E K(n,m)X n Z m p(X n Z m )\ (11) 

n, m— 

For d = 2 we have the most general channel acting on a qubit 

A^ y i (p) = Jr(0, 0)p + tt(1, 0)ApA + tt(1, l)YpY + tt(0, (12) 

For instance, if we take 7r(0, 0) = 1 — p, 7r(l, 0) = 7r(l, 1) = tt(0, 1) = |, the channel , becomes the standard 
symmetric qubit depolarizing channel. 

To justify the choice of Weyl's error model, we point out that the quantum codes employed here are designed to 
error-correct arbitrary quantum errors such as A-errors, Z-errors and combinations of the two (Y-crrors). Since the 
Kraus decomposition of the Weyl channel is defined in terms of powers of these aforementioned error operators, it 
certainly constitutes a natural test-bed where quantifying the performance of the selected qudit codes. Such a test-bed 
turns also out to be very general, with the possibility of encompassing physically relevant scenarios. 

Hereafter, when we consider the two-dimensional Weyl channel, we mean to take into consideration the following 
channel parametrization 

tt (n, m) = f ttx (n) ttz (jn) = ir (n) tt (m) , (13) 



with, 



In this way we have 



tt(1) =tt(-1) d = p, tt(0) = f I -p. (14) 



A^iyi (P) = (1 - p) 2 p + (1 - p) P XpX^ + p 2 Y P Y^ + p(l - p)Z P ZK (15) 

Then, we can also consider the possibility of having asymmetric X and Z errors' probabilities [17H22] . In such a case 
the probabilities tt (n, m) are defined as follows, 

^(asymmetric) ^ m) def ^ ^ ^ ^ ^ (lg) 

with ttx and ttz not identical. In particular, we shall consider 

vrx(l) = K7r z (l), (17) 

with nz (1) = f p and k € [0, 1/p] to guarantee that up is a valid probability value. Notice that for n = 1 we recover 
the symmetric case. 

Finally, when considering block codes on n qubits, the error map simply becomes A^,® 71 acting on p € ©("H®"). 



B. Entanglement Fidelity 

Entanglement fidelity is a reliable performance measure of the efficiency of quantum error correcting codes |23j . 
Suppose a two dimensional code C is such that C C % with dimc'H = N (here 'H can be either 7if n or 'Hd, hence N 
is either 2™ or d). Then consider errors taking place through a CPT map A : &{T~L) — > © (H) (which can be either 
A.( 2 )®" or A^) written, in terms of the Kraus decomposition, as 

k{p) = Y J A kP Al 

k 



() 



To recover the errors by means of the code C a recovery operation must be applied according to the syndrome 
extraction. Suppose it is described by a CPT map 1Z : & (%) —>&(%) 



If the code is effective, we expect that the map resulting from the composition of A and 1Z restricted to the subspace 
C, namely [1Z o A]| c , will be close to the identity map idc on C. In order to evaluate this closeness we can consider a 
state p = trc ("01 written in terms of a purification |?/>) € C ®C and see how well entanglement (between C and the 
reference system identical to C) is preserved by means of 



aide (M W-l) 



This is the entanglement fidelity [T^] for the map [KoA].. 

In terms of the Kraus error operators, T can be rewritten as [21] 



T(p, [K°A] {c )=J2\tel R 3 A k\\c 

3,k 



Finally, choosing a purification described by a maximally entangled unit vector 6 C ® C for the mixed state 
P= dl4c J C' we obtain 



J |c 



(18) 



This is the expression we will use in the following. 



IV. QUDIT CODES FOR WEYL ERRORS 



In this Section, we analyze in details how the qudit codes devised in Sec. II work on the Weyl channel for d = 18 
and d = 50 and determine their performance by means of the entanglement fidelity. 

Our main motivation to use the qudit codes with d — 18 and c? = 50 is that they represent the lowest-dimensional 
perfect qudit systems where a two-dimensional quantum systems (a qubit) can be encoded and protected against 
arbitrary shift errors of the form X n Z m by one and two units, respectively. Furthermore, the physical dimensionality 
of these codes is chosen so to be as much as possible comparable with the physical dimensionality of code-spaces 
characterizing well-known standard stabilizer error correction schemes capable of correcting arbitrary single-qubit 
errors and, possibly, few two-qubits errors. For this reason, the five O |H) and seven-qubit [THl E] quantum stabilizer 
codes seem to be a convenient choice. In particular, recalling that a stabilizer code is perfect if all the eigenvalues of 
the generators constitute valid syndromes for correcting an error, it turns out that both the five and the qudit code 
with d = 18 are perfect and require minimal quantum resources for their task. 

A. The d — 18 qudit code 

Encoding. The encoding operation is characterized by, 

H 2 3 |0) ^ |0 L ) € H 18 and, H 2 3 |1) i— > \1 L ) € H 18 , 
where the codewords |0l) and are defined according to ([5|, ^ as, 

\0 L ) d ^ f i [|0> + |6) + |12)] and, \1 L ) d £ i [| 3 ) + |9) + |X5>] , 
respectively. The stabilizer group S of this code is generated by the two error operators X 6 and Z 6 . Here r\ = r 2 = 3. 
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A simple calculation shows that provided we restrict our focus to the two-dimensional code space c( d - 18 ) c His, 
the following identities hold, 

I = X 6 = A 12 , X = X 7 = X 13 , X 2 = X 8 = X 14 , X 3 ee x- 3 = X 9 = X 15 , 



A 4 = X- 2 = A 10 = X 16 , X* ee x- 1 = A 11 = X 



1 _ vll _ vl7 



and, 



I = Z 6 = Z 12 , Z = Z' = Z u \ Z 2 = Z 6 = Z 14 , Z 3 ee Z~ 6 = Z" = Z 



7 _ r?13 ryl _ ryd, _ rr\i rr 3 _ ry-3 _ ryQ _ rylb 



Z 4 ee Z- 2 = Z 1Q = Z m , Z 5 ee Z- 1 = Z 11 = Z 17 . 

Therefore the total number of error operators to be considered amounts to be 36 rather than 18 2 . Then, the Weyl 
channel (111 may be rewritten as, 



:S5 



K^{p) = Y,A kP Al 



k=0 



where we have relabeled the error operators as follows, 

A Q d ^ y^O)/, Ax d ^ y/^TjZ, •••,^18 = y /n(-l,3)X- 1 Z 3 , 



.A 



del 



.34 



vM3,-2)X 3 Z~ 2 , A 35 ^T(^3jx 3 Z 



In the symmetric case, the probabilities n (n, m) are defined like in Eq.(13), where now 

7T (1) = 7T (-1) 7T(2)=TT (-2) ^p 2 , 7T (3) EE 7T (-3) = p 3 , 7T (0) = 1 - 2 P - 2p 2 - p 3 . 

Correctability. According to (JoJ) with r± = ?*2 = 3 the set of correctable errors -^correctable * s gi ven by the following 
9 errors, 



■^correctable = {4>> A X , A 2 , A 6 , A 7 , As, A 12 , A 13 , A U } , 



(19) 



where, 



A d ^ \AF(0r0)7, A x = V^W^Z, A 2 d ^ y/n(0, -l)Z~\ A, d ^ v^O^A, A 7 d ^ y/^Y)XZ, 



A 8 d ^ 0r(l, -l)XZ-\ A 12 d ^ y/K(-l,0)X-\ A 13 d = l ^(-1, l)X~ l Z, A 14 d ^ ^M, - l)*- 1 ^ 1 . 

Recovery P erators. Following the recipe provided in |25j . it turns out that the two 9-dimensional orthogonal 
subspaces V° L and V 1l of % 18 generated by the action of ^correctable on l^i) an<A are gi yen by, 



V° £ =5pan(|i;2 i > = A= 

I 1 k 1 V^k 



\0 L ),y and V lL = Span\\vl L ) = 



with fc G j( d - ls ) d = {0, 1, 2, 6, 7, 8, 12, 13, 14}. The coefficients denote the errors amplitudes where, for 

instance, ^ d = - 1). Notice that ^K^) = <W<*ij with ft, ft' G Z( d = 18 ) and i, j G {0, 1}. Therefore, it 

follows that V° L © V 11, = His- The recovery superoperator 7\L o {-Rfc} with ft e l( d=18 ) is defined by means of [25], 



In the case under investigation, the entanglement fidelity ( 18 ) reads 



j-^is)^ defers) (I j 8xJiRoA (-=i8)] 1 £ £ |tr([iM. 

7 ( 2 ) i=o 



i;J| C (d=i8) 



feeK d 



where, 



[RkA 



def / (0 L \RkAi\0 L ) (Oi|^A|U) 
iJ l c<d = 18) ~ I (l L \RkAi\0 L ) (l L \R k A\l L ) 



After a simple calculation, J 7 ^ -18 ) (p) becomes, 



jr(rf=i8) (p) = i _ 4p 2 - 2p 3 + 4p 4 



4p 5 +p 6 . 



(20) 



We stress that this error correction scheme is effective as long as the failure probability 1 — J r ( d = 18 ) (p) is strictly 
smaller than the error probability p. This implies that the d = 18 -dimensional qudi t c ode is effective only when 
< P ;$ 0.24. Furthermore, we point out that in its range of effectiveness, J r ( d=18 ) (p) in (20) is a monotonic decreasing 
function of p. 



Asymmetric errors 



By repeating the steps of the previous Subsection using Eqs.(16) and ([17]) we get 

2p 2 - p 3 + k 2 (-2p 2 + 4p 4 + 2p 5 ) + k 3 (-p 3 + 2p 5 + p 6 ) 



J 7 , 



(d=18) 



asymmct 



ric (P) 



1 



Eqs.(20) and (21) become, to the leading order in p with p -C 1 

^ d = 18 )(p) P i 1 



i - V, rjg^Lu: (?) i - 2(1 + ^y. 

It results that the presence of asymmetric Weyl errors with k < 1 increases the performance of the correction scheme. 
This can be understood by noticing that as soon as n — > the noise model becomes classical-like and the errors 
become of a single type, namely Z type, hence more easy to correct (this limiting case is similar to that investigated 
in Ref.[5]). On the contrary, for k > 1 the performance of the code is lowered by errors asymmetries. 



P<1 



(21) 



The d — 50 qudit code 



The encoding operation in this case is characterized by, 

R 2 3 |0) ^ |0 L ) g H 50 and, U 2 3 |1) ^ |1 L ) g H 5(h 
where the codewords \0l) and jl^) are defined according to ([5]), ([6| as, 



|0i 



def 1 

= 7! 



1 10) + |20) + |30) + |40)] and, |lj 



def 1 

= 7! 



+ |15) + |25) + |35) + |45)], 



respectively. The stabilizer group S of this code is generated by the two error operators X 10 and Z 10 . Here ri = r 2 = 5. 

Indeed, it can be shown that provided we restrict our focus to the two-dimensional code space C , ' d=50 - ) C H 50 , the 
set of all (50 2 ) non-normalized errors can be reduced to 100 operators. Following the same line of reasoning of the 
previous Subsection, it is possible to find the recovery superoperator. After some algebraic calculations, we arrive at 
the following expression for the entanglement fidelity 



jr(d=50) (p) = 1 _ 4p 3 _ 4p 4 _ 2p 5 + ^ + 8p 7 + 8p 8 



4p 9 +p 10 . 



(22) 



We emphasize that this error correction scheme is effective as long as the failure probability 1 — J^( d=5 °) (p) is strictly 
smaller than the error probability p. This implies that the d — 50 -dimensional qudit code is effective only when 
< P ^ 0.43. This p-range of effectiveness is larger than that of the d = 18-dimensional qudit code. Furthermore, 
comparing the p-expansions of (20) and (22) to the leading orders for p <C 1, it follows that 



^= 18 ) (p) 



1 - 



4p 2 < 1 - 4p 3 .F( d = 5 °) (p) 



From the above equation, it follows that the d = 50-dimensional qudit code outperforms the d = 18-dimensio nal q udit 
code in the p-range where both error correction schemes are effective as expected. Moreover, J r ( d=5 °) (p) in (22) is a 
monotonic decreasing function of the error probability parameter belonging in its range of effectiveness. 
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Asymmetric errors 



By taking probabilities as in Eqs.(16) and (17) we obtain for the d = 50 qudit code 



r. 



(d=50) 



asymmetric (?) = (l " ^ ~ ^ ~ P 6 ) + + V + V + 2/) 

+ k 4 (-2p 4 + V + 4p 8 + 2p 9 ) + k 5 (-p 5 + 2p 8 + 2p 9 + p w ) 



(23) 



Eqs.(22| and (23) become, to the leading order in p with p <C 1 

p<1 



J-(d=50) (p) 



i-v, CL.(p)^i-2(i+^- 



Also in this case, for k < 1 the presence of asymmetric Weyl errors increases the performance of the correction scheme, 
while for k > 1 the performance of the code is lowered. 



Furthermore, comparing (21) and (23) we get, to the leading order in p with p« 1. 



t: 



(d=18) 



asymmct 



Tic (p) « 1-2(1 + «V. ^ 



-(d=50) 
asymmet 



P<1 



ric (p) 1 - 2(1 + «V 



That is, the d = 50 qudit code outperforms the d = 18 one for any k € [0, 1/p]. 



V. COMPARISON WITH BLOCK (STABILIZER) CODES 

In this Section, for the sake of comparison, we quantify the performance of the standard five-qubit stabilizer code 
[[5, 1, 3]] [HI [5] and the seven-qubit stabilizer CSS (Calderbank-Shor-Steane) code [[7, 1, 3]] [THl HJ on the tensor 
product of Weyl channels (12) by means of the entanglement fidelity [l2l [24] . 



A. The Five-Qubit Stabilizer Code 

Encoding. The [[5, 1, 3]] code is the smallest single-error correcting quantum code. Of all quantum codes that encode 
1 qubit and correct all single-qubit errors, the [[5, 1, 3]] is the most efficient, saturating the quantum Hamming bound. 
It encodes 1 qubit in 5 qubits. The cardinality of its stabilizer group S is |«S| = 2 5_1 = 16 and the 5 — 1 = 4 group 
generators are given by [T3"] . 

{X1Z2Z3X4, X2Z3Z4X5 , X1X3Z4Z5, Z1X2X4Z5} . 

The distance of the code is dc = 3 and therefore the weight of the smallest error AjAk that cannot be detected by 
the code is 3. Finally, we recall that it is a non-degenerate code since the smallest weight for elements of S (other 
than identity) is 4 and therefore it is greater than the distance d = 3. The encoding operation for the [[5, 1, 3]] code 
is characterized by, 

n 2 3 |o) 1 ^ |o L ) e nf 5 and, n 2 3 \i) ^ |i L ) e nf 5 , 

where the codewords |0l) and \1l) are defined as [T3"] . 

|00000) + 1 11000) + |01100) + |00110) + |00011) + |10001) - |01010) - |00101) + 

- |10010) - |01001) - |10100) - |11110) - |01111) - |10111) - |11011) - |11101) 



|o> -> |o L ) d = f \ 



and, 



11)^11,) = 



|11111) + |00111) + |10011) + |11001) + |11100) + |01110) - |10101) - |11010) + 
- |01101> - |10110) - |01011) - 100001) - 1 10000) - |01000) - 100100) - 100010) 



respectively. 

Then, the action of the channel A^ 2 ^ 5 on p £ & (Hf 5 ) can be written as, 



A^(p) = £ A kP At 



k=0 



(24) 
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where, according to the notation of Section II, A k oc X [01020,^,0^0^) Z (&1&2&3&4&5) ■ 

Correctability. The sixteen weight zero and one quantum error operators in ( 24 ) are given by, 

A = y/poh <Z>l2®h®h® I5, M = \fp\~X x ® I 2 <g> h <E> h <8> A 15 = ^/p\^h ® h ® h ® h ® Z 5) 



where the coefficients p\ with I — 0,..., 15 can be easily deduced from the above distribution of errors and Eq.(13l, 

Pq = (1 -P) w i Pi =p(l -P) 9 , •••,P6 =p 2 (l -p) 8 , P15 =p(l -p) 9 - 

It is straightforward, though tedious, to check that, for the above given errors of weight zero and one we have: 

S (^AjA k )j 0, with I, k e {0, 1,..., 15} , 

where S (^A^Akj is the error syndrome of the error operator A\ Ak defined according to ^ as, S (^A^Akj == 

H^ 5 ' 1,3 ^v A i Ak . The quantity H" 5 ' 1 ' 3 " is the parity check matrix for the five-qubit code while v A ^ Aj _ is the vector 

in the 10-dimensional binary vector space F 2 W corresponding to the error operator A\A k . 
Hence, the set of correctable error operators is given by, 



A 



correctable 



{A Q , A t , A 2 , A 3 , A±, A 5 , Ae, A 7 , A s , A g , A 10 , An, A 12 , A 13 , A 14 , A 15 } C A, 
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where the cardinality of A defining the channel in ( 24 ) equals 2 

Recovery Operators. All weight zero and one error operators satisfy the error correction conditions [15] , 

P c A\A k P c cx P C) l,ke {0, 1,..., 15} 
where Pq is the orthogonal projector operator (Pq = Pq and, Pq — Pq) on the code space C^ 5,1 ' 3 ^ defined as, 

Pc = \0l) (0 L | + |U)(U|. 

The two sixteen-dimensional orthogonal subspaces V° L and V 1l of Jif 5 generated by the action of -^correctable 011 
and I It) are given by, 

V° L =Spon||«2*) = ^M°£>}' and V ' L = s P an {\ v l L ) = 

with k = 0, 1,..., 15. Notice that (v t l L \v J l , L ) = S^Sy with I, I' € {0, 1,..., 15} and i, j e {0, 1}. Therefore, it follows 
that V° L © V 1l = "H® 5 . The recovery superoperator 1Z -n- {Ri} with I = 1 ,...,16 is defined by means of [2"5] , 

^ = 10^(^1 + 11^(^1. 



Finally, the composition of this recovery operation TZ with the map A^ 2 )® 5 (p) in (24) yields, 

2 10 -1 16 



(25) 



fe=0 1=1 



Entanglement Fidelity. Here we want to describe the action of IZoA^® 5 in (25) restricted to the code subspace 
C[[5i 1 >3]]_ Note that the recovery operators can be expressed as, 

R l+1 = R x A = (|0 L ) (0 L | + \l L ) A 
vPz VP; 



with i g {0,..., 15}. Recalling that in this case A/ = A], it turns out that, 



(i L \Ri+iA k \j L ) 



{*l\0 l ) {0 L \A\A k \j L 



i L \lL){lL\A}A k \ jL 



We now need to compute the 2x2 matrix representation [ii;Afc]| C of each RiA k with I = 0,..., 15 and k = 0,..., 2 10 — 1 
where, 

w . , def / (0 L |i? /+1 A fc |0 L ) (0 L \R l+1 A k \l L ) 
[U l+1 A k \ lc - 1 (lL | i?/+iAfe | 0i) (i^i^Afclu) 
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For I, k — 0,.., 15, we note that [i?;+i^4fe]| C becomes, 

, fOilAj^lOi 

[-fli+l-AfeJlC 







(l L \A\A k \l L/ 
while for any pair (I, k) with I. = 0,..., 15 and k > 15, it follows that, 



'vAk 



1 
1 



(0 L \Ri +1 A k \0 L ) + (l L \Ri +1 A k \l L ) = 0. 
We conclude that the only matrices [i?;Afc]| C with non- vanishing trace are given by, 



ic 



with s = 1,.., 16. Therefore, the entanglement fidelity (18 1 can be written as 

^[[5,1,3]] (p) ^[[5,1,3]] Q/ 2x2 ,^oAW) =^ 2 E E|tr([^ fe ] |c ) 
V " ' ^' fc=0 Z=l 



and results, 



jr[[5,i.3]] (p) = i _ 40 P 2 + 200/ - 490/ + 728/ - 700/ + 440/ + 175/ + 40/ - Ap 
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(26) 



We remark that this error correction scheme is effective as long as the failure probability 1 — J 7 !! 5 ' 1 * 3 !] (p) is strictly 
smaller than the error probability p. This implies that the five-qubit code is effective only when < p < 2.9 x 10 -2 . 
Finally, we emphasize that this block-encoding scheme is less efficient than the previously-mentioned qudit-encoding 
schemes as it appears in Figure [TJ 



1. Asymmetric errors 



By taking probabilities as in Eqs.(16) and (17) we obtain for the five-qubit code 



[[5,1,3]] 



asymmetric 



(p) 



: 1 - 10/ + 20/ - 15/ + 4/ 
k (-20/ + 80p 3 - 120/ + 80/ - 20/) 
n 2 (-10/ + SOp 3 - 220/ + 280/ - 170/ + 40/ ) 
k 3 (20/ - 120/ + 280/ - 320/ + 180/ - 40p 8 ) 
/ (-15/ + 80/ - 170/ + 180/ - 95/ + 20/) 
k 5 (4/ - 20/ + 40/ - 40/ + 20/ - 4/°) . 



(27) 



We stress that unlike the finding uncovered in [17j . asymmetries in the considered Weyl noisy channel do affect the 
performance of the five-qubit code quantified in terms of the entanglement fidelity. This difference is ultimately a 
consequence of the fact that while in |17j it is assumed error probabilities px, Py and pz all linear in the error 
probability p (although weighted with different coefficients ax, ay and az with ax + ay +az = 1), here we assumed 
Py quadratic in p while keeping both px and pz linear in p. 

From Eqs.(|26|) and (271 we obtain, to the leading order in p with p« 1, 



T 



[[5,1,3]] 



40/ 



T. 



[[5,1,3]] 



asymmetric 



10(1 + k) 2 P 2 



Hence, also in this case, the presence of asymmetric errors increases the performance of the correction scheme for 
k < 1, while for k > 1 the performance of the code is lowered. 
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Hp) 




FIG. 1: The quantities J r(18) (solid thick line), J^ 5 ' 1 - 3 '! (solid thin line), J^ 7 ' 1 * 3 ]' (dash line), are plotted vs. the error 
probability p with < p < 2.6 x 10~ 2 . Within the scale resolution of the graph the curve for J r( - 50 ^ appears to coincide with 
the top horizontal axis. 

B. The Seven-Qubit Stabilizer Code 

The [[7, 1, 3]]-CSS code encodes 1 qubit in 7 qubits. The cardinality of its stabilizer group S is |<S| = 2 7_1 = 64 
and the set of 7 — 1 = 6 group generators is given by [13] , 

1X4X^X^X7, X2X3XQX7 , X1X3X5X7, Z^Z^ZqZ-j, Z2Z^ZqZj, Z1Z3Z5Z7} . 

The distance of the code is dc — 3 and therefore the weight of the smallest error A)Ak that cannot be detected by 
the code is 3. Finally, we recall that it is a non-degenerate code since the smallest weight for elements of S (other 
than identity) is 4 and therefore it is greater than the distance dc = 3. The encoding for the [[7, 1, 3]] code is given 

by \m, 



|o> lew = 



10000000) + 10110011) + 11010101) + 11100110) + 
+ 10001111) + 10111100) + 11011010) + 11101001} 



and, 



|i)-HU) 



(^) 



llllllll) + 11001100) + 10101010) + 10011001) + 
+ 11110000) + 11000011) + 10100101) + 10010110) 



Following the same line of reasoning of the previous Subsection we can compute the correctable errors by means of 
jyll 7 ! 1 . 3 ]]^ the parity check matrix for the seven-qubit code. Finally, after determining the recovery operators, it can 
be shown that the entanglement fidelity reads, 



jr[[74,3]] (p) = l- 42p 2 + 140p 3 + 231/ - 2772/ + 9240p 6 - 18 216p 7 + 24255p 8 - 22 792p 9 



15 246p 10 - 7140p n + 2233p 12 - 420p 13 



36p 
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(28) 



Observe that the seven-qubit code is effective for < p < 2.6 x 10 2 . Comparing the p-expansions of (26) and (28) 
to the leading orders for p <C 1, it follows that 



jr[[5,i,3]] 



1 - 40p 2 > 1 - 42p 2 P « 1 



and, in addition, the p-range of applicability of the five-qubit code is larger than that of the seven-qubit code. Thus, 
for the symmetric Weyl channel considered, the five-qubit code outperforms the seven-qubit code. However, we shall 
see that this ordering does not hold when considering asymmetric scenarios. 
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Asymmetric errors 



In the asymmetric scenario of Eqs.(16) and (171 we obtain for the seven-qubit code 

r[[7,l,3]] 



P. 



... ,„ ic (P) = 1 - -'/'' ! + 70/ - 105/ + SI,/' - 35/ + (i// 

+ k (-21/ + 126/ - 315/ + 420/ - 315/ + 126/ - 2l/) 
+ k 2 (-21/ + 126/ - 315/ + 420/ - 315/ + 126/ - 2l/) 
+ k 3 (-35/ + 420/ - 1785/ + 3850/ - 4725/ + 3360/ - 1295/ + 210p 



10 \ 



k 4 (105/ - 1050/ + 4095/ - 8400/ + 9975/ - 6930/ + 2625/° - 420/ 1 ) 
k 5 (-126/ + 1155/ - 4284/ + 8505/ - 9870/ + 6741p 10 - 2520/ 1 + 399p 
k 6 (70/ - 609/ + 2184/ - 4235/ + 4830/° - 3255/ 1 + 1204/ 2 - 189/ 3 ) 



-15/ + 126/ - 441/ + 840/° - 945/ 1 + 630/ 2 - 23l/ 3 + 36p 



14\ 



(29) 



From Eqs.(28) and (29) it follows, to the leading order in p with p -C 1, 

^ (/ 1 - 42/, (p) *g 1 - 21(1 + .)/. 



Once again, for k < 1 it results that the presence of asymmetric errors increases the performance of the correction 
scheme, while for k > 1 the performance of the code is lowered. 



Furthermore, by comparing (27) and (29) it follows, to the leading order in p with p <C 1, 

«Lc (P) 1 - 10(1 + «) V > 1 - 21(1 + Ktf ^ (p) , 1-1 

[[5,1,3]] 



asymmetric 



(p) 



i - io(i + «) v < i - 21(1 + «y j-ISU w , - 



> i.i 



Thus, with respect to the noise model discussed in Ref.|17j. we conclude that here the comparison between five-qubit 
code and seven-qubit code is slightly more involved. 



VI. FINAL REMARKS 



In this article, we discussed how to protect a qubit embedded into a qudit from both amplitude and phase errors 
occurring in the discrete phase space. A code has been devised using stabilizer formalism and its performances 
compared with those of common block codes for a general Weyl noisy quantum channel allowing symmetric and 
asymmetric error probabilities. 

Specifically we have considered the d = 18 and d = 50 qudit stabilizer codes together with five and the CSS 
(Calderbank-Steane-Shor) seven-qubit quantum stabilizer codes. The performances of these codes were quantified by 
means of the entanglement fidelity as function of the error probability. 

We uncovered that qudit codes have an enormously wider (by approximately an order of magnitude) range of 
applicability in the error probability. Furthermore, already the d — 18 qudit code outperforms the five and seven- 
qubit block codes for symmetric errors (see Fig[I]). Our theoretical analysis leads to the conclusion that the qudit 
codes with d = 18 and d = 50 outperform the common five and CSS seven-qubit stabilizer codes. This in principle 
allows one to save space resources (since d = 18 < dime Hf 5 = 32 << dime Hf 7 = 128), however one should also 
account for the difficulties in implementing qudit systems, an issue that seems to be nontrivial and not quite settled 
yet. For an overview of the experiments performed for producing quantum optical qudits, we refer to [26] . 

The performance of qudit code are also robust against asymmetries in errors' probabilities. In fact, restricting our 
analysis to n > 1.1, it results that the d = 18 qudit code outperforms the seven-qubit code until strong asymmetries 
come into play, as can be seen by comparing (21 ) with (29) to the leading order in p with p <C 1, 



[[7,1,3]] 



asymmetric 



(p) 



p<i 



1 - 21(1 + n)p 2 < 1 - 2(1 + k 2 )/ 



P«l j~{d=l 



asymmetric 



(p) , K < 



21 + V593 



11.34 



•FlL!±l trlc (p) P ^ 1 - 21(1 + n)p 2 > 1 - 2(1 + « 2 )/ P ^ ic (P) , - > « H-34 



asymmctr 



Comparative results for the various codes performances in presence of asymmetries are graphically represented in 
Figure [2] 
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error probability p for k — 4 < 11.34 (left) and for k = 20 > 11.34 (right). Within the scale resolution of the graphs the curves 
for -T^fymmotric appear to coincide with the top horizontal axes. 



The different uncovered behaviors in the four error correcting schemes employed in this article can be ascribed to 
the fact that the errors in V n ®s (or V H »7) are fundamentally different from those in V-Uig an d 7 > Hw 

Finally, it could be interesting to consider the presence of correlations between X and Z errors in the qudit code. 
These can be introduced as follow 

7T (n, TO,) = f (1 — /i) 7T (n) 7T (to) + (J,8 ni m 7T (?Tl) , 

where fi with < ji < 1 represent the degree of correlation. Following the very same line of reasoning provided in 
Section IV, it can be shown that for instance the entanglement fidelity becomes, 

Ho7t W) (p) = (! - V ~ 2 P 3 + 4 P 4 + V + P 6 ) + /•* (2p 2 + p 3 - 4/ - 4p 5 - p 6 ) . 
It then results (dT^ 18) /dfi) P « 2p 2 > 0, that is memory effects lead to better performances. The reason 

V / p— const. 

is that in the limit of very strong correlations fi —> 1, only one type of error (namely Y = XZ) takes place. As such, 
this case shows similarities with the case of asymmetric errors with k — > 0. 

In conclusion, we are strongly motivated by our investigation to believe that encoding a qubit into a qudit can be a 
useful approach in quantum coding. We are aware of the difficulties in realizing and controlling qudit systems even of 
low dimensionality, however we have witnessed a lot of progress along this direction recently. Quantum optical qudits 
can be generated by means of experimental schemes based upon interferometric set-ups, orbital angular momentum 
entanglement and, biphoton polarization [26] . For instance, in the interferometric scheme employed in 27 , high 
symmetry and maximally entangled qutrits are realized with a fidelity up to 0.985 as the superposition state of the 
three possible paths of a single photon in a three-arms interferometer. Therefore, the realization of the discussed 
qudit codes seems not futureless. 
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